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Abstract. We define an approximate version of the Fourier transform on 2^ 
elements, which is computationaUy attractive in a certain setting, and which may 
find apphcation to the problem of factoring integers with a quantum computer as is 
currently under investigation by Peter Shor. [SHO] 

Fourier Transform 

Notation: Let L be a positive integer. Let a, c be L-bit integers. The binary 

representations of a, c are 



a 



!^Q^aj2', c — E, 




Define the L-bit integer h as the reversal of c. 





In terms of the binary representations. 




Whenever j -\- k > L, we have u;(^'^'^ = 1, so that we can drop those terms from 
consideration: 



{FFT) 




Notice that all computations are in the field Q(a;(2i)). 
Hadamard Transform 

The Hadamard transform looks like a Fourier transform defined over . It suits 
my purposes, pedagogically, to reverse the indexing on the output of the Hadamard 
transform, and get the transform 



(Normally the exponent would be "EjajCj, but we reverse the indexing to bring out 
the similarity with the ordinary FFT.) Rewrite this as 

= -l=E„X„exp (^^E}-^ajCL_,_j2^-'^ 

1 /^27ri 

— ^^^^SaA^exp ^-^Ijo<j,k<L-l;j+k=L-l(ljCk2 

noting that exp(27ri2^~^/2^) = —1. Since the sum is restricted to those values of j, k 
satisfying j + k — L — 1, we can replace 2^~^ by 2-^+'^ and obtain 

1 /27ri 
(HT) Yc = -j=EaXaexp (^^So<j- fc<L-i;i+fc=i,-iajCfc2^+'' 

Approximate Fourier Transform 

Comparing the two formulas (FFT) and (HT), wc find that the only difi^erence is 
in the limits on j + /c : in (FFT) the range isO<j-|-/i;<L — 1, while in (HT) the 
range isL — l<j-|-A;<L — 1. 

This leads us to define an Approximate Fourier Transform ( AFFT) , parameterized 
by an integer m : 

1 /27ri \ 

(AFFTjn) = ^^^/=So-'^aexp |^-^Eo<j,fe<L-l;L-m<j+fe<L-l%Cfe2'' J 

When m = 1 this is the Hadamard transform (suitably indexed); when m = L it 
becomes the ordinary Fourier transform. 

Since j + k > L — m, the argument of "exp" is some multiple of 27ri2^~"*/2^ = 
27ri/2"^, so that AFFT is defined over Q(c<J(2"i)). 

The argument of "exp" in AFFT differs from that of FFT by 



The magnitude of this difference is bounded by 

2^ 

If L = 500 and m = 20, this bound is about 3/1000. So the matrix entries of AFFT 
differ from those of FFT by a multiphcative factor of exp(ie) where |e| < 27rL2~™ = 
3/1000. Thus if AFFT is used in place of FFT in Shor's factoring work [SHO], it leads 
to an overall error of a fraction of a degree in each phase angle, and less than one 
percent decrease in the magnitude of the probability of each desirable final state. 
Calculating the AFFT 

Start with the description of the Fast Fourier Transform as taken from [KNU, page 
291, section 4.3.3]. I have replaced A, t, s, k by X. a. 6, L, respectively, and numbered 
the passes from L — 1 down to 0, to correspond to the bit being manipulated. 

* Initialization. Let X[^l(aL_i, Oq) = Xa, where a — (aL_i...ao)2 (the binary 
representation). 

* Pass J, J = L - 1, L - 2, 1, 0. (Numbered downwards!) Set 

* X^-^\bL-i, bL-2, bj, aj_i, ao) := 

X[-^+i](6i_i, bj+i, 0, aj_i, ao) + 
^(6,6,+,...6,_,o...o). X xi'+'KbL-i,...,bj+ul,aj_„...,ao) 

We wish to compute the FFT quantum mechanically. At the outset, X^^^(aL-i, oo) 
represents the amplitude of the state where L electrons have spins ol-i, Oq, respec- 
tively, with "1" representing "up" and "0'' representing "down". Each succeeding 
X^-^^{bL-i, ao) represents the amplitude of the state of these same L electrons. The 
transform is performed by a sequence of two-electron interactions. 

On Pass J, multiply the amplitudes X['^+^1(6l_i, 1, oj-i, Oo) (with a 1 in 

position J) by the phase shift ^;(o^J+i---''i-iO---0)2, This correspond to the following two- 
bit operations. For each K, J+1<K<L — 1, use an interaction between electrons 
J and K to multiply the amplitude of those states with a 1 in both positions J and 
K by the factor 

Call this transformation Qjk- 

Then apply the unitary transformation 

71(1 -1) 



3 



to the electron J. Call this transformation Pj. So for L 
K is K = 2, and we have 



3, J = 1, the only value of 







" 1 



































1 



































1 


o;^ 































1 


1 

































7/9 

V z 











1 





cu^ 





= 
























1 





cu^ 



























1 





a;6 

































1 





a;6 _ 














1 


1 


















" 1 























1 


1 


















1 

















1 


-1 
























1 

















1 


-1 













X 











1 




















1 


1 
















1 




















1 







1 
















1 

















1 


-1 
























a;2 














1 







-1 

























1 

72 



In general one would have L — 1 — J of the two-bit interactions Qjk on pass J, 
corresponding to different values of K. The entire 3-spin FFT is depicted in the 
Appendix. 

So the FFT matrix is expressed as a product of unitary matrices. For example 
FFT on 4 electrons is 

-foQoi Qo2 Qo3 -Pi Q 12 Q 13 -P2Q23-f3 ■ 

If there are L electrons then there are L matrices Pj and L{L — l)/2 matrices Qjk- 
For our approximate AFFT, we simply delete those matrices Qjk with K > 

J + m. So the AFFT is again unitary, and easily computed with one-bit and two-bit 

operators. It requires about Lm two-bit operations. 
Quantum computation 

Shor [SHO, page 12] suggests first developing a state 

1 



^a=0 



a > 



where q ~ is a product of small prime powers, which will enable him to do a 
mixed-radix FFT later. By contrast, we suggest setting q = 2^ ^ 5n^. Second, he 



4 



computes x"'{mod n), where x,n are integers computed classically, so that the state 
becomes ^ 

— pS^^gla, x"' > 

Then he does the Fourier transform, sending a to c with amplitude exp{2Triac/q). 
This leaves the machine in state 

^S;^~Loexp(27rmc/g)|c,a;" > . 

We see that the radix-2^ Fourier transform is directly implementable as 2-spin 
interactions, as opposed to the operations required by Shor. 

We can improve still further, by doing our approximate Fourier transform instead 
of the Fourier transform. Notice that on Pass J of AFFT computation, we use 
interactions between bits J and K,J<K<J + m. So bit i^T of the output index, 
bx, does not participate in any interaction after pass J — K — m. (Remember we are 
numbering backwards, so pass K — m is m passes later than pass K.) Similarly, bit 
K of the input index, a^^, does not enter into the computation until pass J — K. 

So we propose rearranging the computation in the following way. 

* Start with 7/ = 1 in an L-bit quantum register 

where you will compute x". 

* For each J = L - 1, L - 2, 2, 1, : 

* Place the electron J in state 

^(|0>+|l» 

corresponding to the two possible values of oj. 

* Compute 

y :— y{x^ )'^°'-'\mod n) 

reversibly, in the register allocated for y. 

* For K = J + \, J + 2, min( J + m — 1, L — 1), apply operation Qjk- 

* Apply operation Pj. 

* If J < L — m, measure bit = CL-j-m from the output 

of pass J of the AFFT computation. 
(It will not enter any more interactions.) 

* End (For each J = L - 1, L - 2, 2, 1, 0). 

* Measure the remaining bits 6^-2, &o- 
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* End algorithm 

A possible advantage of this arrangement is that the electron in position K need 
only maintain coherence for m passes of the computation, although the rest of the 
system still has to maintain coherence for a longer time, so this advantage might be 
less than it appears at first blush. 

A definite advantage is in the computational complexity. Shor's proposal, using 
a mixed-radix Fourier Transform with q k, 5n^ the product of small prime powers, 
appears to require about (logn)^ elementary operations (spin-spin interactions). The 
radix-2^ FFT requires only (logn)^ elementary operations. The AFFT requires only 
(log n) (log log n + log 1/e) operations, where a final precision of e is required. So the 
Fourier transform is no longer the bottleneck of the computation. 

Parallel implementation 

Several steps of the AFFT can be parallelized in the quantum implementation; 
this might further speed up the computation time, and increase the likelihood of the 
state remaining coherent until the computation is done. 

We use, for an example, the FFT on 5 electrons, with operations proceeding right 
to left: 

FFT = P0Q01Q02Q03Q0AP1Q12Q13QUP2Q23Q24P3Q34P4 

We can interchange the order of any two operations which do not involve any of the 
same electrons; alternatively, we can do such operations in parallel. At time step 
K — 8,7, 1, 0, let us perform Pi ii I + I = K , and Qu ii I + J = K. Steps that 
are performed in parallel are displayed within square brackets, vertically aligned, and 
again proceeding right to left: 



FFT = Po<9i 



01 



Pi 
Q02 



Q12 

Qo3 



P2 
Ql3 
Qo4 



Q23 
Ql4 



P3 

Q24 



Q34P4 



We used 9 time steps here; for an L-electron system we will use 2L — 1 time steps. 
This parallel implementation looks a lot like "systolic arrays," [MC, chapter 8, 
section 8.3], and suggests directions for physical implementation. 
References 

[KNU] Donald E. Knuth, volume 2. The Art of Computer Programming, Volume 
2: Seminumerical Algorithms. (Addison- Wesley, Reading, MA, 2nd ed., 1981) 

[MC] Carver Mead and Lynn Conway, Introduction to VLSI Systems. (Addison- 
Wesley Reading, MA, 1980) 



6 



[SHO] Peter W. Shor, "Algorithms for Quantum Computation: Discrete Log and 
Factoring," manuscript, 1994. Proceedings of FOCS 1994. 
Appendix 

We write out in full the FFT on 3 electrons. Note that the rows are numbered in 
bit-reversed order (04261537), corresponding to the index of b rather than c. 

FFT = P0Q01Q02P1Q12P2 — 
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UJ OJ UJ UJ UJ UJ UJ UJ 



LO U) U) LO LO LO 
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LO 



UJ UJ 

4 2 
UJ UJ 
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Uj' 
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